Improving fidelity of continuous-variable teleportation via local operations 
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We study the Braunstein-Kimble setup for teleportation of quantum state of a single mode of 
optical field. We assume that the sender and receiver share a two-mode Gaussian state and we 
identify optimum local Gaussian operations that maximize the teleportation fidelity. We consider 
fidelity of teleportation of pure Gaussian states and we also introduce fidelity of the teleportation 
transformation. We show on an explicit example that in some cases the optimum local operation is 
not a simple unitary symplectic transformation but some more general completely positive map. 
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I. INTRODUCTION 

Quantum state teleportation is undoubtedly one of the 
most exciting developments in the rapidly growing field 
of Quantum Information Processing. In quantum tele- 
portation, the information about the teleported quan- 
tum state is transferred from the sender, Alice, to the 
receiver, Bob, via dual classical and quantum EPR chan- 
nels The latter is established via an entangled state 
shared by Alice and Bob. The teleportation protocol goes 
as follows: Alice carries out a Bell-type measurement on 
the state she wants to teleport and her part of the shared 
entangled state. She sends the result of her measurement 
via classical channel to Bob, who applies to his part of 
entangled state a transformation which depends on the 
classical information received from Alice. 

The teleportation is perfect and Bob recovers an exact 
copy of the state teleported to him by Alice only if the 
quantum channel is ideal maximally entangled state. If 
we deal with qubits represented by polarization states 
of photons, then we can employ pair of polarization- 
entangled photons generated by means of spontaneous 
parametric down-conversion, where the entanglement is 
almost perfect [|J 0. However, in case of continuous 
quantum variables ||, ^) , an ideal EPR channel is an un- 
physical infinitely squeezed state. In quantum optics, the 
available resource is a two-mode squeezed vacuum state 
with some finite degree of squeezing [|[ ||. Moreover, the 
parts of the entangled state must be distributed among 
Alice and Bob, e.g., through optical fibers. This trans- 
mission inevitably introduces losses and noise and the 
entangled state shared by Alice and Bob will be some 
mixed state in general. 

An important question is whether one can somehow 
improve the quality of the teleportation by means of local 
operations on the parts of the shared entangled state. Re- 
cently, this problem has been studied for teleportation of 
qubits and it was shown that local transformations may 
indeed be helpful [j?], |j| |j| . Moreover, it was demonstrated 
that the optimum local transformation that maximizes 
the average teleportation fidelity need not be simple uni- 
tary transformation, but some completely positive (CP) 
ma P PI- I 11 other words, it may be advantageous to 
let the parts of the shared quantum state interact with 



local ancillas. 

In this paper, we investigate how to improve the fi- 
delity of teleportation of continuous quantum variables 
by means of local operations on sender's and receiver's 
side. The first steps in this direction were already taken. 
Bowen et al. showed that in certain cases the fidelity of 
teleportation of coherent or squeezed states may be im- 
proved when Alice and Bob locally apply squeezing trans- 
formations to their parts of the shared quantum state 
Jlpj l . Kim and Lee considered an asymmetric mixed quan- 
tum channel and showed that in that case the fidelity of 
teleportation of coherent states may be enhanced when a 
local transformation accompanied by decohcrcncc is ap- 
plied to one part of the quantum channel . 

To make the problem tractable, we restrict ourselves to 
the class of trace-preserving Gaussian CP maps H, [l3| . 
These maps preserve the Gaussian shape of the Wigner 
function of the transformed state. The restriction to 
Gaussian CP maps is very reasonable from the experi- 
mental point of view, because these maps can be imple- 
mented in the laboratory as a unitary symplectic trans- 
formation (linear canonical transformation of quadrature 
operators) on the signal mode and auxiliary modes ini- 
tially prepared in some Gaussian states. In quantum 
optical setups this can be done with the help of phase 
shifters, beam splitters and squeezers. Gaussian CP 
maps were recently applied to description of cloning of 
continuous quantum variables [fl3[ . Another recent paper 
discussed the conditions under which a given two-mode 
shared Gaussian state can be transformed into another 
Gaussian state by means of local Gaussian CP maps . 

The paper is organized as follows. In Sec. II we 
will briefly describe the Braunstein-Kimble teleportation 
setup and we will derive compact formulas for fidelities of 
teleportation of any pure Gaussian state. We shall also 
introduce a fidelity for the teleportation operation itself. 
It will turn out that this latter fidelity can be interpreted 
as a fidelity of entanglement swapping. In Sec. Ill we will 
briefly review the properties of Gaussian CP maps and 
we will derive optimum local Gaussian CP map which 
maximizes a chosen teleportation fidelity. We shall con- 
sider two scenarios: in the first case the transformation 
is applied only on one side, in the second case both Al- 
ice and Bob may locally apply some CP maps. In Sec. 
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IV we present an example of our optimization procedure. 
Finally, Sec. V contains conclusions. 



II. FIDELITIES 

We shall consider the Braunstein-Kimble setup for 
teleportation of a single mode of optical field ||. The 
quantum channel between Alice and Bob is established 
via two-mode entangled state pab fully described by 
its Wigner function Wab(xa,Pa, %b,Pb)- Alice mixes 
the mode whose state she wants to teleport with her 
part of entangled state on balanced beam splitter and 
she carries out a homodyne detection on each output 
mode thereby measuring two commuting quadratures 
X+ = (x in + x A )V2 and P_ = [p in - p A )/V2. After 
receiving the measured values of X + and P_ from Al- 
ice, Bob displaces his part of entangled state as follows: 
xb — > xb + \/2X+, ps — > Pb + V2P-. We assume ideal 
homodyne detectors on Alice's side and a zero coherent 
component of the entangled state p A B (mean values of all 
quadratures xa,b, Pa,b vanish). Under these conditions 
the resulting state on Bob's side possesses the same co- 
herent component as the original state teleported to him 
by Alice and the teleportation is invariant under displace- 
ment transformation. 

Ide et al. showed that the fidelity of continuous 
variable (CV) teleportation can be improved by optimiz- 
ing the gain g in the modulation of the output field whose 
quadratures are displaced by the amount gX+ and gP- . 
However, in this case the teleportation is not in general 
invariant under displacement transformation and, for in- 
stance, the fidelity of teleportation of coherent state |a) 
depends on the intensity |a| 2 . Here we keep the gain g 
fixed and improve the teleportation fidelity by suitable 
local transformations of the shared entangled state. As- 
suming fixed gain g = y/2, the relation between input 
and output Wigner functions of the teleported state is 
given by convolution Jl6|] 

/oo 
K(x 2 -X\ , p 2 -pi ) W in [xi , pi )dx 1 dp 1 . 
-OO 

(1) 

In order to express the kernel K it is convenient to rewrite 
Wab as a function of the variables x± = xa ± xb and 
P± = Pa±Pb, 

W A b(xa,Pa,x b ,Pb) =Wab(x + ,p + ,x~,p-). (2) 

With the help of Wab we can write 

1 f°° 

K{x + ,p_) = - \ W A B(x+,P+,x~,~P-)dx_dp + . 
J — 00 

In what follows we shall assume that the shared quantum 
state pab is two-mode Gaussian state. This is reasonable 
assumption since this class of states can be prepared in 
the lab. It is computationally convenient to deal with 



characteristic function of this state, defined as Fourier 
transform of the Wigner function, 



1 f°° 

Wab{t) = ——7 / wab (<z) exp(iq • r) d 4 q, 



(4) 



where r = (x a ,Pa,x b ,Pb) and q = (£a,Va,£b,Vb) are 
real vectors. For Gaussian state with vanishing coherent 
component we have |17|| 



wab {q) = exp 



1 ^ 7 
--qT A Bq 



(5) 



The elements of the covariance matrix Tab are given by 

F ABtij = (AnArj) + (ArjAn), (6) 

where Arj = rj — (rj) (note that if the coherent compo- 
nent of the state vanishes then (rj) — 0). We can express 
the real covariance matrix Fab hi terms of three 2x2 
matrices A, B, and C, 



ab 



A C 
C T B 



(7) 



Here A and B are covariance matrices of the single modes 
on Alice's and Bob's side, respectively, and C contains 
the inter-modal correlations. 

Now consider teleportation of a pure single-mode 
Gaussian state with covariance matrix D. Since the 
teleportation is invariant under displacement transforma- 
tion, all states with the same covariance matrix but dif- 
ferent coherent components are teleported with the same 
fidelity. It thus suffices to consider state with vanishing 
coherent amplitude, whose characteristic function reads 



Win(qii 



exp 



1 ^ T 

-q iu Dq in 



(8) 



where q m = (£i n ,??in)- Fidelity of teleportation of a pure 
state can be calculated as an overlap integral of input 
and output Wigner functions over the whole phase space 

/oo 
W in (x,p)W out (x,p)dxdp. (9) 
-00 

After making use of the formulas (Q) and (||) , expressing 
all Wigner functions as Fourier transforms of the charac- 
teristic functions and carrying out all integrals, we arrive 
at a compact formula for the teleportation fidelity, 



F 



(10) 



Vdet E ' 
where the matrix E reads 

E = 2D + RAR T + RC + C T R T + B (11) 

and 



R 



1 

-1 



(12) 
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Besides the fidelity of teleportation of certain class of 
Gaussian states, one can introduce the fidelity of the tele- 
portation process itself. How this can be accomplished 
becomes clear when one notices that the teleportation 
transformation ([!]) is a trace-preserving CP map [p7s[ . 
Any CP map can be represented by positive semidcfinitc 
operator x on a Hilbert space which is a tensor product 
of the Hilbert space of input states Tt and Hilbert space 
of output states JC This representation is not only 

mathematical, the state x can be actually prepared in 
the lab if we first prepare a maximally entangled state 
on Hilbert space H® 2 and then apply the CP map to one 
part of the entangled state. In case of CV teleportation, 
the maximally entangled state is the EPR state 

Wepr = ^-S( Xl - x 2 )S( Pl + p a ) (13) 

and the teleportation of one part of that state can be 
interpreted as an entanglement swapping |2(], ^l| . Hence 
the fidelity we obtain in this way is the fidelity of entan- 
glement swapping of the EPR state. 

Formally, the CP map that transforms input density 
matrix p m onto output density matrix p ut can be written 
as a partial trace over the input Hilbert space, 



Pout = Tr H [xp[ n <8> tic}- 



(14) 



In our case, the Wigner function Wt e \ of the teleportation 
CP map xtoi is closely related to the kernel K because 
the convolution ([[]) is essentially the partial trace jl^ ) 
rewritten in terms of Wigner functions, 



Wtd = — K(x 2 -Xi,p2 +pi) 
Zn 



(15) 



Notice the change of sign in front of pi which reflects the 
transposition in Eq. (|l4|). The ideal teleportation is an 
identity map represented by the EPR state (|l3|). Now 
since the CP maps are represented by positive semidefi- 
nite operators and since the ideal transformation is rep- 
resented by a pure state (|l~3|), we can calculate the fidelity 
between the ideal and actual teleportation as fidelity of 
these two states 22 . Thus we can write 



/oo 
W EPR (xi, Pl ,X 2 ,P2) 
-oc 

xWtei{xi,pi,X2, P 2)dx 1 dpidx2dp 2 - 

(16) 

On inserting the explicit formulas (|l^) and (|l^) into Eq. 
( |l6| ) we obtain 



/oo 
dxdp. 
-OO 



(17) 



We can see that T x is infinite, as could have been ex- 
pected since we work in infinite dimensional Hilbert 
space. Nevertheless, the fidelity ( |l7| ) can be renormal- 
ized. If we drop an infinite constant proportional to Dirac 
delta function and multiply by 2n, then we obtain 



If we insert the explicit formula (||) for kernel K into Eq. 
(|l8|), then we find that 

/oo 
Wab(x,p, —x,p)dxdp. (19) 
-oo 

For Gaussian quantum channels (^) , this formula sim- 
plifies to 



T = 



Vdet E' 
where the matrix E' reads 

E' = RAR T + RC + C T R T + B. 



(20) 



(21) 



The expression for the fidelity T is a rather special case 
of the formula for the fidelity of teleportation of pure 
Gaussian states ( |l0| ) where we set the covariance matrix 
D equal to zero. Of course, this means that T is un- 
bounded. Nevertheless, T is a good measure of the qual- 
ity of teleportation. For instance it can be shown that 
T > 1 only if the state pab is entangled (see Appendix). 
In particular, if pab is two-mode squeezed vacuum state 
parametrized by squeezing constant r then one gets 



T — exp(2r), 



(22) 



hence the fidelity monotonically exponentially grows with 
the squeezing. 



III. OPTIMUM LOCAL GAUSSIAN CP MAP 

Our task is to maximize the fidelity of teleportation (ei- 
ther F or T) by means of local Gaussian trace-preserving 
CP maps on Alice's and Bob's side. We shall consider 
two scenarios: in the first, simpler scenario the CP map 
is applied only on Bob's side while in the second case 
both Alice and Bob may locally apply some CP maps. 

Gaussian CP maps are those maps for which the 
Wigner function of the corresponding operator x has a 
Gaussian form. The teleportation with Gaussian quan- 
tum channel is an example of Gaussian CP map p8| . 
The Wigner function representing single-mode trace- 
preserving Gaussian CP map reads 



27r 2 A/delG 



exp (-AT- T G- x Ar) , (23) 



where S and G are real 2x2 matrices, moreover, G is 
symmetric positive semidefinite matrix, 



Ar 



Srt 



(24) 



and 



( Zout \ 


< = ! 


Zin \ 


V Pout / ' 




v -Pin j 



T = 2ttX(0,0). 



(18) 



are column vectors of output and input quadratures, re- 
spectively. Since we deal with Gaussian states whose 
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form is invariant under Gaussian CP maps, it suffices to 
provide rule for transformation of the covariance matrix 
I\ The relation beteween input and output single-mode 
covariance matrices Ti n and r out is given by a simple 
linear map [O 



sr in s 2 



G. 



(25) 



The map ( p3| ) is completely positive if and only if S 
and G satisfy an inequality |lq] 



G + iS-iS , S5 T > 0, 



where 



1 

-1 



(26) 



(27) 



The condition ( |26| ) can be derived as follows: the Wigner 
function ( jga ) must represent a positive semidefinite oper- 
ator, which imposes constraint on the covariance matrix 
G p3|. Namely, the matrix 



M 4 . 



Gij + iAr^Arj], 



(28) 



must be positive semidefinite, where [, ] stands for com- 
mutator. Making use of canonical commutation relations 



for the quadratures x lni p m and Xq-^u 



one arrives af- 



ter some algebra at the inequality (|26j). 

Assume now that a Gaussian CP map ( p3| ) is applied to 
Bob's part of shared two-mode state pab- This modifies 
the covariance matrix Tab, 



ab 



A CS T 
SC T SBS T + G 



(29) 



The maximization of the fidelity then amounts to the 
minimization of the determinant 

V = dot [2D + RAR T + RCS T + SC T R T + SBS T + G] 
_ (30) 
under the constraints ( Pq ) . Recall that on setting D = 
we obtain as a special case the fidelity of entanglement 
swapping j2C|). 

We divide the optimization of the CP map into two 
steps. In the first step we find optimum G for a given 
matrix S and then we shall optimize over all possible 
matrices S. Since the matrices S and G have altogether 
seven independent elements, 



G 



( .9ii 


.912 \ 


s = 


Sll 


Sl2 \ 


V 512 


.922 ) ' 




\ S21 


S22 J 



(31) 



we have to find a global minimum of a function of seven 
real variables under the constraint (p6[), which can be 
equivalently expressed as 



and 



.9ii > 0, 322 > 



3ii322 - 3i2 - (1 - s ) 2 > 0, 



(32) 
(33) 



where s = S11S22 — si2S2i- We introduce a short-hand 
notation for the elements of matrix 



2D + RAR 1 + RCS 1 + SC T R T + SBS 



a 7 
7 P 



(34) 

Notice that this matrix is, by definition, positive semidef- 
inite, and its elements a, (3, 7 are functions of Sij. Thus 
we can write the determinant (pCfl) in a compact form, 



(35) 



V = (a + gu)(P + 322) - (312 + l) 2 



It is always optimal to choose "extremal" matrix G 
that satisfies the inequality ( P3|) as an equality. Indeed, 
if a sharp inequality holds in ( |33|) , then we can reduce the 
value of diagonal elements 311 and 322 until the equality 
is reached in ( |33"| ) and this would obviously reduce also 
the value of T>. Hence we can write 



312 



= ±^311322 - (1 - s) 2 



(36) 



and insert into (|35|). Furthermore, we can see that it is 
optimal to choose the sign of 312 the same as the sign of 
7 and we have, 



V = (a+.9ii)(/3+322Mv / 3ii.922 - (1 - s? + \l\f - (37) 
Upon solving the set of two nonlinear extremal equations 

(38) 



dV 



dgi 



0. 



dV 



dg2 



= 0, 



we find that the optimum matrix G is proportional to 
the matrix (34), 



G 



V 



a 7 
7 



(39) 



On inserting the elements of the optimum G back into 
Eq. (|37]) we finally obtain 



V= f|l-*| + Va/3-7 2 



(40) 



Now T> is a function of four variables sn, S22, S12 and 
S21 and we have to find its global minimum. In general, 
such optimization is a hard task and can be solved only 
numerically. However, we shall see that when making 
some assumptions we will be able to solve this problem 
analytically. 

It is well known that by means of local symplectic 
transformations it is possible to bring any two-mode co- 
variance matrix Tab into tridiagonal form (24), 



AB 



a a 

a c 2 

a b 

c 2 b 



(41) 



It suffices to consider Gaussian quantum channels for 
which all the matrices A, B and C in (|7]) are diagonal. 
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Further assume that also the covariance matrix D of the 
teleported state is diagonal, D = diag(d n , d 22 ). Note 
that this assumption is not a serious restriction, because 
any D can be diagonalized by means of reversible sym- 
plectic transformation. In this case it can be shown that 
the necessary conditions on extremum 



dV 

dsi2 



0. 



dV 

ds 2 i 



0. 



(42) 



are satisfied when s\2 = s 2 \ = 0. One has to be a bit 
careful here because there is an absolute value in Eq. (Efl) 
and three cases must be distinguished: (i) 1 — s > 0, (Ii) 
1 — s < 0, and (iii) s = 1. In all cases, the conditons on 
extremum ( |42"| ) are satisfied when S12 = S21 = 0. 

We are thus lead to make the hypothesis that if the 
matrices A, B, C, D are diagonal, then the optimum G 
and S are also diagonal. When S is diagonal, then 7 = 
0, and the matrix elements a and f3 become quadratic 
functions of su and S22, respectively, 



a(sn) = 2dn + a + 2c\S\\ + &sf 1; 
P(s22) = 2c?22 + a - 2c 2 s 2 2 + bs^- 



(43) 



For the sake of notational simplicity we define x = sn 
and y = S22 and we must minimize the function 



f{x,y) = |1 - xy\ + y/a(x)P(y). 



(44) 



The extremal equations are obtained by setting the par- 
tial derivatives of f(x,y) equal to zero, 



x = ±4 



/ a(x) 



(by -<%), y = ±\ 



a(x) 



(bx+a), (45) 



where the signs + and — correspond to the cases when 
1 — xy > and 1— xy < 0, respectively. From the product 
of the formulas for x and y, we can express y in terms of 
x, 



C2(bx + a) 
V ~ x{b 2 -l) + bci 

Substituting this formula back into the second Eq. 
and squaring that equation, we arrive at 



(46) 



c\ol(x) = [x(b 2 - l) + b Cl ] 2 fi 



C2(bx + Cl) 

x{b 2 - 1) + bci 



(47) 



This is a quadratic equation for x and can be solved an- 
alytically. In this way we identify all potential minima 
outside the boundary xy = 1. It remains to localize 
minima on the boundary where y = 1/x and we must 
minimize the function a(x)(3(l/x). The condition on ex- 
tremum 



d_ 

dx 



a(x)(3 ( - 
1 x 



= 



(48) 



reduces to quartic equation for x. Upon solving this equa- 
tion we get positions of all possible minima on the bound- 
ary, i.e., we determine all potentially optimum symplectic 
transformations. 

Let us now consider a more general protocol, where 
both Alice and Bob are allowed to apply some local Gaus- 
sian CP maps. To make the problem tractable, we do 
not assume any communication between Alice and Bob 
at this stage, hence they both apply their local opera- 
tions independently. Furthermore, we shall assume that 
all relevant matrices are diagonal, hence we shall seek the 
optimum two-mode CP map in the form 



(u 

v 

0x0 

V y , 



(49) 



G = 



9A.ll 














9A,22 














5B.11 














9b, 22 



(50) 



The covariance matrix Tab transforms according to 



AB 



ST 



AB* 



G. 



(51) 



From Eq. (|33|) where the equality should hold and where 
.912 = 0, we obtain the following relations between the 
elements of the optimum matrix G, 

3a,ii3a,22 = (1-iw) 2 , 9b.ii9b.22 = i^-xy) 2 . (52) 
The determinant T> can be expressed as 

(1 — uv) 2 (1 — xy) 2 ~ 



V = (a + g A .n + 9b, 11) + 



9A.11 



9b,u 



(53) 




This function attains its global minimum when 
fot 

9A.ii = |1 - uv\J —, gB,ii = \l-xy 

On inserting these expressions back into Eq. (|5 
V = (|1 - xy\ + |1 -uv\ + V"/3) 2 

where a and (3 are functions of four real variables 
u,v,x,y, the elements of matrix S. In general, the mini- 
mum of the function (|55| ) must be found numerically. In 
what follows we shall focus on the fidelity of entangle- 
ment swapping and we shall see that in this case one can 
find the global minimum analytically. The square root 
of the determinant ( |55| ) that we must minimize reads in 
this case (du = c?22 = 0) 

f(u,v,x,y) = |1 - uv\ + |1- xy\ 

+ [(u 2 a + 2uxcx + x 2 b)(v 2 a - 2vyc 2 + y 2 b)] 1/2 . 

(56) 
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This is actually a function of only three variables. This 
becomes apparent when we make the following substitu- 
tions 



0.65 



TO -> W, X — > x/v, 

The function (|5q) then reads 



y -> yv- 



(57) 



f(w,x,y) = |1 - w| + |1- a;y| 

+ [(u) 2 a + 2wa;ci + x 2 fo)(a - 2yc 2 + y 2 &)] 1/2 . 

(58) 



After another substitution x = qw the function (|58|) be- 
comes a linear function of w: 



f(w, q, y) = |1 - w\ + |1 - wgy| 

+ [(a + 2gci + g 2 6)(a - 2yc 2 



1/2 



(59) 



From the linearity of (|59j) it is clear that the extrema are 
localized at points, where one absolute value is equal to 
zero. Hence we have to consider three different possibil- 
ities: 

(i) w = 1, no operation is applied on Alice's side and 
a CP map is applied on Bob's side. 

(ii) wqy = 1, a symplectic transformation is applied 
on Bob's side. However, this symplectic transformation 
can be in our case "absorbed" into CP map on Alice's 
side, hence another possibly optimum strategy is to do 
nothing on Bob's side and to apply a CP map on Alice's 
side. 

(hi) w = 0, this means that both Alice and Bob throw 
away their parts of shared quantum state and replace 
them with vacuum states. Clearly, this strategy is opti- 
mum if the quantum channel is not in entangled state, 
because with vacuum state at both sides one gets max- 
imum fidelity T obtainable without the aid of entangle- 
ment, J~ max, class — 1- 

One may object that the substitution x — wq is prob- 
lematic when w — and x 0. However, a detailed anal- 
ysis reveals that if one of the four parameters u,v,x,y is 
set equal to zero and the three remaining parameters are 
optimized, then we once again arrive at the above listed 
alternatives (i)-(iii). 

The strategies (i) and (ii) represent a CP map on only 
one side, while nothing is performed on the other side. It 
was shown above that these optimum one-sided CP maps 
can be found analytically. The search for optimum CP 
map would thus consist of three parts: find one-sided 
optimum Gaussian CP maps on Alice's side, on Bob's 
side and also consider replacement of the shared quan- 
tum state with vacuum state and choose the optimum 
alternative that yields maximum fidelity. 



IV. EXAMPLE OF OPTIMIZATION 

To illustrate how the optimization works in practice, 
let us assume that the covariance matrix Tab has the 



F 0.55 




0.45 



FIG. 1: Fidelity of teleportation of coherent state for bo = 0.5 
and variable squeezing r. The solid line shows the maximum 
fidelity achievable via local CP map on Bob's side and the 
dashed line shows the maximum fidelity achievable via local 
symplectic transformations on Bob's side. Both curves coin- 
cide when the squeezing is higher than the threshold rth- 



tridiagonal structure given by Eq. 
elements read 



ill) and the nonzero 



a = 1 + 2 sinh 2 r 
6=1 + 2 sinh 2 r + b Q 



Ci = — sinh(2r), 
C2 = sinh(2r). 



(60) 



For bo = we recover the covariance matrix of pure two- 
mode squeezed vacuum state and the quantum channel 
is in a mixed state for any b > 0. 

Let us analyse how the fidelity of teleportation of co- 
herent state can be improved by means of local trans- 
formations on Bob's side. Since e?n = d<zi = 1 and also 
ci = -C2 [c.f. Eq. @], the solution of Eq. @ sim- 
plifies considerably because the functions a and /3 are 
identical. The optimum x and y are equal, x — y, and 
the two roots of Eq. (|47|) read 



C2 



f 



X2 



C-2 



1 



(61) 



Furthermore, the quartic equation ([48]) for optimum sym- 
plectic transformation splits into two quadratic quations 
that have only two real roots x = ±1. 

We must evaluate the fidelity of teleportation of co- 
herent state for all these potentially optimum transfor- 
mations on Bob's side, and choose the maximum value. 
The resulting fidelity is plotted in Fig. 1 for bo = 1/2 
and a variable degree of squeezing r. ft turns out 
that if the squeezing is lower than certain threshold 
r t h = — [ln(l — £>o)]/2, then the optimum transformation 
on Bob's side is a CP map with x given by x± in Eq. 
d6l|). The parameter x\ grows from zero for r — to the 
value x\ = 1 that is attained when r — r t h- For higher 
squeezing, the best strategy is to do nothing, i.e., the 
optimum operation is a symplectic transformation with 
x = y = 1. The optimum CP map for r < r t ^ is a 
simple damping process which can be implemented with 
the help of a beam splitter with amplitude transmittance 
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FIG. 2: The same as Fig. 1 but the fidelity of entanglement 
swapping T is plotted. 



t = C2/Q) — 1) whose two input ports are fed with Bob's 
part of entangled state and a vacuum state, respectively. 
This transformation reduces the noise represented by 60 
in Eq. (|60|), which in turn improves the teleportation 
fidelity. 

With the help of CP map on Bob's side, the fidelity 
of teleportation of coherent state is always larger than 
the maximum fidelity 1/2 achievable without the aid of 
entanglement. On the other hand, if we allow only for 
unitary symplectic transformations on Bob's part of the 
state, then there is a region of squeezing where the max- 
imum achievable fidelity is lower than 1/2. This example 
clearly illustrates that in certain cases it is advantageous 
to couple the shared state to the local environment J?], |) . 

Similar results are obtained for the fidelity of entan- 
glement swapping T . In this case we can optimize over 
all local Gaussian CP maps on both Alice's and Bob's 
side, because the problem reduces to the optimization of 
a Gaussian CP map on only one side, as discussed in the 
previous Section. For our specific example, the optimum 
CP map is actually the same as that for the fidelity of 
teleportation of coherent state. Also the dependence of 
the fidelity T on r is qualitatively similar to that shown 
in Fig. 1, see Fig. 2. In particular, there is a region 
where !F > 1 if Bob applies the optimum CP map, but a 
restriction to local symplectic transformation results in 
T < 1. 



CONCLUSIONS 



We have restricted ourselves to the class of local trace- 
preserving Gaussian completely positive maps and we 
have shown that in this case the optimization problem 
can be solved analytically. We have demonstrated on a 
simple example that the optimum local operation need 
not be a unitary symplectic transformation but some 
more general CP map. 
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APPENDIX A 

Here we prove that the inequality T < 1 holds if the 
shared quantum state Wab{xa,Pa,xb,Pb) is separable. 
Our starting point is the the formula jis|): 



/oo 
Wab(x,p, —x,p)dxdp. 
-00 



(Al) 



Density matrix pab of any separable state can be written 
as a convex mixture of product states, 



Pab 



E 



P, Paj <s Pi:..- 



(A2) 



where pj > and J2j Pj = 1 • Formula ( |A2j ) implies that 

W A b(xa,Pa,x b ,Pb) = / IPjWAj (xa,Pa)Wb,j (xb,Pb)- 

(A3) 

If Wbj(xb,Pb) is Wigner function of the quantum state 
Pbj, then Wbj(—xb,Pb) is a Wigner function of the 
quantum state p because the transformation x — ► — x 
and p — > p is the transposition. For separable state QA2D, 
the formula (Al) thus reduces to 



F = ^2Pj Tr iPA,jpB,j}, 



(A4) 



In this paper we have shown that one can improve 
the fidelity of teleportation of continuous quantum vari- 
ables by means of local operations on the sender's and 
receiver's parts of the shared entangled state pab (quan- 
tum channel). We have considered the fidelity of tele- 
portation of pure Gaussian states and we have also in- 
troduced a fidelity measure for the teleportation transfor- 
mation. The latter fidelity was interpreted as a fidelity of 
entanglement swapping of infinitely squeezed EPR state. 



where we used that 



/oo 
W A {x,p)W B {x,p)dxdp. (A5) 
-00 

The Schwarz inequality implies that Tr[pA,jpB j] < 1 and 
with the help of normalization of pj we finally obtain 



T<1. 



(A6) 
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The fidelity T = \ forms a boundary between classical in- illustrates the essential and central role of the entangle- 
formation transfer and quantum teleportation. The state ment in the teleportation. 
Pab must be entangled in order to achieve T > 1. This 
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